Symplectic Schemes for Birkhoffian system * 



Hongling Su Mengzhao Qin 

CAST (World Laboratory), 
Institute of Computational Mathematics 
and Scientific/Engineering Computing, 
\ Academy of Mathematics and System Sciences, 

O ■ Chinese Academy of Sciences, Beijing, 100080, China. 



O 

(N 

C 

Hi 

m 



O 
O 



February 7, 2008 



Abstract 



A universal symplectic structure for a Newtonian system including 
nonconservative cases can be constructed in the framework of Birkhof- 
fian generalization of Hamiltonian mechanics. In this paper the sym- 
plectic geometry structure of Birkhoffian system is discussed, then the 
■ symplecticity of Birkhoffian phase flow is presented. Based on these 

CO \ properties we give a way to construct symplectic schemes for Birkhof- 

fian systems by using the generating function method. 

£3 ; 1 Introduction 

Birkhoffian representation is the generalization of Hamiltonian representa- 
tion, which can be applied to hadron physics, statistical mechanics, space 
^ ■ mechanics, engineering, biophysics and so on^|5]. All conservative or non- 

conservative, self-adjoint or non-self-adjoint, unconstrained or nonholonomic 
constrained systems always admit a representation of BirkhofF's equations^ 
13]. In last 20 years, many researchers have studied Birkhoffian mechanics 
and obtained a series of results in integral theory, stability of motion, inverse 
problem, algebraic and geometric description, and so on. 

BirkhofF's equations are more complex than Hamilton's equations, so are 
the studies of their computational method. In the past, there are not any re- 
sult in the computational methods for Birkhoffian system. The known differ- 
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ence methods are not generally applicable to Birkhoffian system. As the dif- 
ference scheme to solve Hamiltonian system should be Hamiltonian scheme, 
the difference scheme to simulate Birkhoffian system should be Birkhoffian 
scheme, however, the conventional difference schemes such as Euler center 
scheme, leap-frog scheme and so on are not Birkhoffian scheme, so a method 
of how to systematically construct a Birkhoffian scheme is necessary. It is 
the main context in this paper. 

The systems both Birkhoffian and Hamiltonian mentioned in this paper 
are usually finite dimensional situation py, in fact, the definition of infinite 
dimensional Birkhoffian system has not proposed before. Ref. |2j described 
the Algebraic and Geometric profiles of the finite dimensional Birkhoffian 
systems in local coordinates, and general nonautonomous Hamiltonian sys- 
tem is considered as autonomous Birkhoffian system. Ref. jlj developed 
systematically the symplectic schemes for the standard Hamiltonian system 
and general Hamiltonian system on the Poisson manifold which belong to- 
gether the autonomous and semi-autonomous Birkhoffian system. So in this 
paper, we just discuss the nonautonomous Birkhoffian system in detail. 

In section 2, we sketch out Birkhoffian system via variational self-adjoin- 
tness, this sketch shows the relationship between Birkhoffian system and 
Hamiltonian system more essentially and directly, then we give the basic 
geometrical properties of Birkhoffian system. 

Section 3 extends the definitions of i^(z)-Lagrangian submanifold to the 
one with a parameter t such as K(z, i)-Lagrangian submanifold, then we dis- 
cuss the relationship between symplectic mappings and gradient mappings. 

Section 4 constructs the generating functions for the phase flow of the 
Birkhoffian system and gives the method to simulate Birkhoffian systems 
by symplectic schemes of any order. The last section shows an illustrating 
example. A first-order and a second-order schemes for the system of a linear 
damped oscillator are given. 

The Einstein's summation convention is used in the next sections. 

2 Birkhoffian System 

The generalization of Hamilton's equations we shall study is given by 




dB(z,t) dF(z,t) 
dzi ' dt 




(1) 
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following the terminology suggested by Santilli(1978) .2], we called it as 
Birkh- 

off's equation or Birkhoffian system. The function B(z,t) is called the 
Birkhoffian, because of certain physical difference with Hamiltonian which 
is indicated in Ref. 2 . Fi, i = 1,2,- • -,2n, are Birkhoffian functions. 
A representation of Newton's equations via Birkhoff 's equation is called a 
Birkhoffian representation. 

Definition 2.1. Birkhoff 's equations 0j are called autonomous when the 
Fi and B functions are free of the time variable, in which case the equations 
are of the simple form 

They are called semi- autonomous when the Fi functions do not depend ex- 
plicitly on time, in which case we have the more general form 

Birkhoff 's equations are called nonautonomous when both the Fi and B func- 
tions have explicit dependence on time, in which case we rewrite 

dzj dB(z,t) dFi(z,t) 
dt dzi dt 

Here they all have 



Kij (z,t)^- - - ^£1L = 0. (4) 



Kl3 ~ d Zl dzi 

They are called regular when the functional determinant is not null in the 
region considered: 

det(Kij)(U) + 0, (6) 

otherwise, degenerate. 

Given an arbitrary analytic and regular first-order system 

K ij {z,t)-^ + D i {z,t) = Q, i = l,2,...2n, (7) 

From the point of the inverse variational problemj^j, this system is self- 
adjoint iff. it satisfies the following condition in i.e. 

Kij + Kji = 0, 

OKjj dKjk dK ki _ Q 

dz k dzi dzj ' (8) 

dKij dD t dDj . . 

—^T- = -^ hJ,k = 1,2, ...2n. 

dt dzj dzi 
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We now simply introduce the geometric significance of the condition of vari- 
ational self-adjointness[Hl El QUI Ell Q2j- Here the region considered is a 
star-shaped region of points of R x T*M, T*M the cotangent space of 
the M, M a 2n-dimensional manifold. 

Consider first the case for which Kij = Kij{z). Given a symplectic 
structure written as the 2-form in local coordinates 



one of the fundamental properties of symplectic form © is that dfi = 0. 
The geometric significance of the condition of self-adjointness (JBJ is then 
the integrability conditions for 2-form Q to be an exact symplectic form 
coincident with the first two formulas of condition (JHJ). Because the exact 
character of two-form Q implies following structure 



this geometric property is fully characterized by the first two equations of 
condition (JSJ), and we can say that the two- form (|10[) describes the geomet- 
rical structure of the autonomous case (J2J) of the Birkhoff 's equations, even 
it sketches out the geometric structure of the semi-autonomous case. 

For the case of that = Kij(z,t), the full set of condition Q must 
be considered and the corresponding geometric structure can be better ex- 
pressed by transition of the symplectic geometry on the cotangent bundle 
T*M with local coordinates Zi to the contact geometry on the manifold 
R x T*M with local coordinates Zi, i = 0, 1,2, • • 2n, zq = iM. In this case 
more general formulations of an exact contact 2-form persist, although it is 
now referred to (2n+l)-dimensional space, 

2n 



2n 

17 = Kij(z, t)dzi A dz. 
»J=1 



K — —K 



(9) 




(10) 





where 




(12) 



if the contact form is also of the exact type 





(13) 
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the geometric meaning of the condition of the self-adjointness is then the 

integrability condition for the exact contact structure (|13j) . Here B can be 
calculated from 

-P = A + ^ (14) 

dzi dt v ' 

for 

^( A + i>r) = &-( D ' + "arj- (15) 

All the above discussion can be expressed via the following property. 

Proposition 2.1. (Self—Adjointness of Birkhoffian System). Necessary a- 
nd sufficient condition for a general nonautonomous first- order system given 
as above to be self-adjoint in of points of R x T*R 2n is that it is of the 
Birkhoffian type, i.e., 

„ / s dzi ^ , x ,dFj dFi.dzi , x dF(z,t). 

A', J(2 ,*)^ + = - ^ - (V-B(M) + -^). (16) 

The functions F and -B can be calculated according to the rules [3] 

l 

F i = ±Jz r K ji (\z,t)d\, (17) 
o 

i 

r bf- 

B= I Zi .(D i + -^)(\z,t)d\. (18) 
o 

As well as Hamiltonian system, Birkhoffian system provides a symbio- 
sis among variational principle, Lie's algebra and symplectic geometry. |131 
n HI EI Obviously both standard Hamilton's equations and general Hamil- 
ton's equations on Poisson manifold are recovered from Birkhoff 's equations 
as in the particular cases of autonomous and semi-autonomous Birkhoffian 
representations. 

Due to the self-adjointness of Birkhoff 's equations, the phase flow of the 
system (fH)|) conserves the symplecticity, then we get 

= ^(Kijdzi A dzj) = 0. (19) 

It means that if we denote the phase flow of the equations (|16|) with (z, t) , 
then 

Kij(z, t)d% A dzj = Kij(z, t)dzi A dzj, (20) 
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or 

— K(z,t)— = K(z,t). (21) 
oz oz 

In the next sections we will construct the algorithm preserving this geometric 
property of the phase flow in discrete space. 

3 Generating Functions for K(z, £)-Symplectic Ma 
ppings 

In this section we consider general k(z, t)-symplectic mappings and their 
relationship with the gradient mappings and their generating functions. 

Definition 3.1. Denote 

j2n= (-l 0)' j4 - = (-L t)' 

J (Jin \ ~ (K(z,t) \ 

J4 " = V0 -jj' K(z,z,t,t )=[ - K{zM )- 

A 2n- dimensional submanifold L C R An 

L=|^ * e R 4n \z = z(x,t ),z = z(x, t),x eU C R 2n , open set 

(23) 

is a J4 

n - or «^4n* or K(z, z,t,to)-Lagrangian submanifold if 

(T x L) T J 4n (T x L) = (24) 

or 

(T x L) T J 4n (T x L) = (25) 

or 

(T x L) T K(z,z,t,t )(T x L) = 0, (26) 



?2n 



where T X L is the tangent space to L at x. 

A mapping with parameters t and to is z — ► z = g(z,t,to) : R 2 
R 2n which is called a canonical map or a gradient map or a K(z,t)-symplectic 
map if its graph 

r fl = {( * ) eR 4n \z = g(z,t,t ),z = zeR 2n } (27) 

is a J4 

n~ or Jiw or K(z, z,t,to)-Lagrangian submanifold. 
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Definition 3.2. A differentiable mapping g : M — > M is K{z,t)-symplectic, 
if 

^ K(g(z,t,t ),t)^=K(z,t ). (28) 
oz oz 

A difference scheme approximating the Birkhoff's system Mb}) 

z k+1 =g k (z k ,t k + T,t k ), k^O, (29) 

either explicit scheme or worked out from a inexplicit scheme, is called a 
K-symplectic scheme, when g k is K-symplectic for every k ^ 0, i.e. 

|VV W )| = ^,^). (30) 

The graph of the phase flow of the Birkhoffian system Q is g t (z,to) = 
g(z,t,to) which is a K(z, z, t, io)-Lagrangian submanifold for 

gl{z, t ) T K(g\z, t ),t)gl(z, t ) = K{z, to). (31) 

Similarly the graph of the phase flow of standard Hamiltonian system is a 
J4 n -Lagrangian submanifold. 

Define nonlinear transformation with two parameters t and to from R in 
to itself, 



a (Mo) : 



z \ / w \ = ( ai(z,z,t,to) 

z ) \ w J \ a 2 (z,z,t,t ) 



«- X (Mo) : ( t )^( I ) = ( if^Y^ )• (32) 



Denote 



w J \ z J \ a (w,w,t,to) 
a*(z,z,t,t ) 





B a 


C a 


D a 


(A a 


B- 


[c a 


D 



a, 1 (w,ro,Mo)=La ™ • (33) 



a* is the Jacobian of a. Let a be a diffeomorphism from R 4n to itself, 
then it follows that a carries every K-Lagrangian submanifold into a J^ n - 
Lagrangian submanifold, iff. 



J in a* = K, (34) 
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i.e., 

'A a B a \ T f,J 2n \(A a B a \ = (K(z,t) \ 

C a D a J \0 -J 2 n)\C a D a J V -K(z,t )J- { ' 

Conversely a~ 1 carries every J4 n -Lagrangian submanifold into a iT-Lagrang- 
ian submanifold. 

Proposition 3.1. Let M £ R 2nx2n j a be defined as above. Define a frac- 
tional transformation 



a a :M — > M 

M ^N = a a (M) = (A a M + B a )(C a M + D a )~ l 

under the transversality condition 



(36) 



(38) 



\C a M + D a \^0. (37) 

Then the following four conditions are equivalent mutually: 

\C a M + D a \^0, 
\MC a - A a \ 
\C a N + D a \ ^ 0, 
\NC a -A a \ ^0. 

The proof is direct and simple, so it is omitted here. 

Theorem 3.1. Let a be defined as above. Let z — > z = g(z,t,to) be a 
K(z,t)-symplectic mapping in some neighborhood 7Z of R 2n with Jacobian 
g z (z,t,to) = M(z,t,to). If M satisfies the transversality condition in 1Z 

\C a {g(z, t, t ), z, t, t )M{z, t, t ) + D a (g{z, t, t ), z, t, t )\ + 0, (39) 

then there exists uniquely in 1Z a gradient mapping w — ► w = f(w,t,to) 
with Jacobian f w (w,t,to) = N(w,t,to) and a scalar function- generating 
function-cj)(w,t,to) such that 

f(w, t, t ) = 4> w (w, t, t ), (40) 
ai(g(z,t,to),z,t,t ) = f(a 2 {g(z,t,to),z,t,to),t,t ) 

(41) 

= 4>w(a2{g(z, t, t ), z, t, t ),t, t ), 
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identically in z and t, 



N = (A a M + B a ){C a M + D a )~ , (42) 

M = (A a N + B a )(C a N + D a )- 1 . (43) 
Proof. Under the transformation a, the image of the graph T g is 

a(T ^ = {{ Z ) ei?4n l^ = «i(s(^Mo),^Mo), 

w = a2(s(z,Mo),z,Mo) } • 

By the inequality (|3ll|) . 

dw da>2 dz 8ct2 „ ,^ „ ,„,. s 
-K - = -K^r • p — h -p— = C a M + Z> a ^ 0, 45 
az oz oz oz 

so w = ot2(g(z,t,to),z,t,to) is invertible, the inverse function is denoted by 
z = z(w, t, to). Set 

w = f(w,t,t ) = ai(0(2,Mo),z,Mo)U=z(t O ,t,t o ), (46) 

then 

N = JT = + ^TP) = ^ M + B °W° M + D-r 1 - (47) 

au; oz oz oz ow 

Notice that the tangent space to a(T g ) at z is 



dz 



we have a conclusion that a (r g ) is a J4 n -Lagrangian submanifold for 
T z (a(T g )) T J An T z (a(r g )) 

= ((A a M + B a ) T ,(C a M + D a f)j 4n ( cImXdI ) 

= (M T ,/) a Jj 4 „ a ^ M ^ (49) 

= f )=0. 



So 



(A a M + 5 a ) T (C a M + D a ) - (C a M + D a ) T (A a M + B a ) = 0, (50) 



i.e., N = (A a M + B a )(C a M + D a )~ l is symmetric. It implies that w = 
f{w,t,to) is a gradient mapping. By the Poincare lemma, there is a scalar 
function (f>{w,t,to) such that 

f(w, t, t ) = 4> w {w, t, t ). (51) 

The equation ([41 1) follows from the construction of f(w, t, to) and z(w, t, to)- 
Since z(w, t, to)oa2(g(z, t, to), z, t, to) = z, so substituting w = a2(g(z, t, to), 
z,t,to) in the equations (|4*H|) and ([5Tjl . we get the equation (|41j) . 

Proposition 3.2. f(w,t,to) obtained in Theorem lff.il is a/so t/ie solution 
of the following implicit equation 

a l {f{w, t, t ), w, t, t ) = g(a 2 (f{w, t, t ), w, t, t ),t, t ). (52) 

The proof is similar to the above proof. 

Theorem 3.2. Let a be defined as in Theorem lff.il Let w — > w = 

f(w,t,to) be a gradient mapping in some neighborhood 1Z of R 2n with Jaco- 
bian f w (w, t, to) = N(w, t, to)- If in TZ, N satisfies the condition 

\C a {f{w, t, t ),w, t, t )N(w, t, t ) + D a (f(w, t, t ),w, t, t )\ + 0, (53) 

then there exits uniquely in TZ, a K(z,t)-symplectic mapping z — > z = 
g(z,t,to) with Jacobian g(z,t,to) = M(z,t,to) such that 

c^{f(w, t, t ),w, t, t ) = g{a 2 (f(w, t, t ),w, t, t ),t, t ), 

M = (A a N + B a )(C a N + LP)- 1 , (54) 
N = (A a M + B a )(C a M + D a y\ 

Similarly to Proposition I ff. S\ g(z,t,to) is the solution of the implicit equation 

ai{g(z, t, to),z, t, t ) = f{a 2 {g{z, t, to), z, t, t ), z, t, t ). (55) 
The proof is similar to that of Theorem 13. II and is omitted here. 

4 Symplectic Difference Schemes for Birkhoff 's Eq 
uations 

In Section 2 it is indicated that for a general Birkhoff's system, there is 
a common property that its phase flow is symplectic. Through the result 
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in the proceeded section, we construct symplectic schemes for Birkhoff's 
system by approximating the generating functions. 

The Birkhoff's phase flow is denoted by g t {z,to), it is a one-parameter 
group of K(z, i)-symplectic mappings at least local in z and t, i.e., 

g to = identity, g tl+t2 = g h o g t2 , (56) 

here z is taken as an initial value when t = to, and z(z,t, to) = g t (z,to) = 
g(t;z,to) is the solution of the Birkhofhan system 



Theorem 4.1. Let a be defined as above. Let z — > z = g t {z,to) be the phase 
flow of the Birkhoff's system Mb]) . M(t;z,to) = g z (t;z,to) is its Jacobian. 
At some initial point z, i.e., t = to, z = z, if 

\C a (z,z,to,t ) + D a (z,z,to,t )\ ^ 0, (57) 

then there exists, for sufficiently small \t — to I an d in some neighborhood o- 
f z G R 2n , a gradient mapping w — » w = f(w, t, to) with Jacobian f w (w, t, to) 
= N(w,t,to) symmetric and a scalar function- generating function-cp(w,t, to) 
such that 

f{w, t, t ) = w {w, t, t ), (58) 



d 

—(j) w (w, t, t ) = A((j) w (w, t, t ),w, 4>ww(w, t, t ),t, t ), (59) 
at 

dw - dw 

A{w, w, Q^,t, to) = A(z(w, w, t, t ),z(w, w, t, t ), g^,t, t ), (60) 



7 dw d ^ dw d 

ai(g(t; z, to), z, t, t ) = f(a 2 (g(t; z, t ),z, t, t ),t, t ) 
= 4> w (a 2 (g(t;z,to),z,t,to),t,to), 

identically in z and t, 

N = a a (M) = (A a M + B a )(C a M + D a )-\ 
M = a a -i = (A a N + B a )(C a N + D 01 )' 1 . 



(61) 



(62) 



(63) 
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Proof. M(t;z,to) is differentiable with respect to z and t. Condition 
(|57j) guarantees that for sufficiently small \t — to\ and in some neighborhood 
z of z £ R 2n , there is 

\C a (z,z,t,to)M(t;z,t ) + D a (z,z,t,t )\ / 0. (64) 

Additionally the Birkhoffian phase flow is a symplectic mapping, therefore 
by Theorem 13. 1 1 there exists a time-dependent gradient map w = f(w,t,to) 
and there is a scalar function <p(w,t,to), such that 

f(w, t, t ) = <p w (w, t, t ). (65) 



df(w,t,t ) d(f> w (w,t,t ) 

— df— = at (66) 

Notice that z = g(t;z,to) is the solution of the following initial-value 
problem 



f = K -\%t){s7B + % )(z,t) 



z\t=t Q — z 
therefore from the equations 



so 



dw dw dw dw 
~dt ~ ~dt ~dt~dt 

,a dw^. T ^_ u _ OF. da\ dw da 2 



(67) 



w = ai(z, z, t,t ), w = a 2 (z,z,t,t ), (68) 

it follows that 

dw did dz 9 

~dt = 8?' dt + 8f Ql( "' M,j) 

+ + (69) 
^- 1(v S + f) + ^, 



(70) 



Since |^ ^ 0, so w = w{w,t) exists and it is solvable, But it can not be 
solved directly from the transformation a and a -1 . Set 

dw dw 

^''to'^-W (71) 
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then we denote 



did 



dw 



A(w, w, Q^it, to) = A(z(w, w, t, t ),z(w, w, t, t ), g^,t, t ), (72) 
so from the equation (|66|) . we get 



The above theorem is given for the nonautonomous Birkhofhan system, 
which is geometrically different from the autonomous or semi-autonomous 
Birkhofhan system and we explained this point in Section 2. According to 
Ref. [1] , we can easily construct any order symplectic difference schemes for 
the later two. Because of the simplicity of the ordinary geometry structure in 
autonomous or semi-autonomous case, the transformation a we need in these 
cases is free of parameter t, and the corresponding Birkhofhan is completely 
Hamiltonian, so 




(73) 



Q.E.D.. 



dw 
~dt 



dw 
~dt 



did dw 
~dt^iJ 




B w (z(w,w),t)) 



(74) 



resonantly we have the Hamilton-Jacobi equation 




(75) 



or for semi-autonomous cases there is 




(76) 
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We describe the k-th. order total derivative of A with respect to t as 

oo 00 

D k t A = d^A(E(t ~ t,) l ^ +l) ) + d^A(£(t - Wfi^H 

i=0 i=0 

00 00 
dtd^ w A(J2(t ~ to)¥ w k ~ 1+l) + dtd^A^it ~ *>)V&T 1+i) )+ 

i=0 i=0 
k k—m k~m—n 

^2 c r^2 c k- m ^2 ^w^ww^t 1 ^ (77) 

m=0 n=l 1=1 hi-i hhn + 

ji-i \-ji=k-m 

00 00 

- t )Vi hl+<) , ■ ■ ■, X> ~ t°Y4>t +l) , 

i=0 i=0 

00 00 

^ ( t-t o) V^),---,^(t-toM^ ) ), 

i=0 i=0 

then at the point of t = to, the total derivative of A is as 

D k t A{^\wA { l,hM = d (j>w A4>^ ) + d (Pww A^i+ 
dtd.A^-V + d t d, w A^ 1] + 



k k—m k—m—n 



n pi 1 , a™ 

(78) 



^2 C ™ ^-"^ X] d 4>w d 4>ww d i 

m=0 n=l 1=1 h 1 H hhn+ 

jiH \-ji=k-m 

A(^\w,At ,t )^tK-- ;^\&K-- 

Theorem 4.2. Lei .4 and a be analytic. Then the generating function 
<t>wa,A( w ,t,to) = <f>w(w,t,to) can be expanded as a convergent power series 
in t for sufficiently small \t — to\ 

00 

4> w {w,t,t ) = 5>-t )Vl fc Wo), (79) 

k=0 

(k) 

and (jr w , k ^ 0, can 6e recursively determined by the following equations 

</><£Hw,to) = f(w,t ,t ), (80) 

$Hw^ = A(<j>$\w,<l>W,t ,t ), (81) 



b k w +1 (w,t ) = j^-^D k A(^\w,^lt ,t ). (82) 
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Proof. Differentiating the equation (f?9|) with respect to w and t, we 

get 

oo 

4> ww (w,t,t ) = Y,(t-to) k ^wl(w,t ), (83) 



fc=0 

and 



ft 

-g f Mw, t, t ) = J2(k + m~ t ) k ^ +1) (w, t ). (84) 
By the equation (j58|) . 

4>Z(w,to) = 4> w (w,to,t ) = f(w,to,t ). (85) 

Substituting the equation (|83|) in A(w, w, ^^,t,to), and expanding A in 
t = to, we get 

A(<p w ,w, (f) ww ,t, t ) = A{f{w, t , t ),w, f w (w, t ,to),to,t )+ 

f:^(t-to) k D^\ w ,^i,to,t ). (86) 

k=l 

Using the equation (|59f) and comparing the equation (|5B]) with the equation 
l|Hij>. we get the equations (jSU and (JH2J. 

For the autonomous and semi-autonomous the difference is that A is re- 
placed with Bikhoffian B, which makes the expanding of generating function 
4> much more easy. 

Through Theorem (|4.1|) and ()4.2I) . the relationship between the Birkhof- 
fian phase flow and the generating function <p(w,t,to) is established. With 
the aid of this result, we can directly construct k{z, t)-symplectic difference 
schemes. 

Theorem 4.3. Let a be given as Theorem \J^. i| ), and A analytic. For 
sufficiently small r > as the time-step. Take 

m 

^(w,to + r,to) = ^2T i ^(w,to), m = l,2,..., (87) 

i=0 

where 4>w are determined by equations Then ^^(w,^ +t, to) 

defines a K(z,t)-symplectic difference scheme z = z k — > z k+1 = z, 

ai(z k+1 ,z k ,t k+1 ,t k ) = ^i m) (a 2 (z k+1 ,z k , t k+1 ,t k ),t k+1 ,t k ) 

of m-th order of accuracy. 
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Proof. Because of that 



\C a (z,z,t ,t ) + D a (z,z,t ,t )\ ^0, (89) 

so by Proposition 13.11 \C a N + D a \ ^ 0, where N = cj> ww (wo,to,to) = 
^ww(wo, to, to), and wq = a(z, z, to, to). Thus for sufficiently small r and 
in some neighborhood of wq, there exists 

\C a N {m \w,to + T,to)+D a \ ^0, (90) 

where 

(w, to + r,t ) = ^ (w, to + T, to). (91) 

By Theorem l3.21 ipw (w, to+r, to) defines a K(z, i)-symplectic mapping 
which is expressed in the equation ([55 j). Therefore the equation ([%%]) deter- 
mines a m-th order K(z, i)-symplectic difference scheme for the BirkhofHan 
system (fl6|) . 

5 Example 

In this section we give an example to illustrate how to get a symplectic- 
structure-preserving scheme for a non-conservative system which can be ex- 
pressed in Birkhoffian representation. 

Consider the system of a linear damped oscillator 

f + vr + r = 0. (92) 

We introduce a gradient function p satisfying p = r, then a Birkhoffian 
representation of this system is given by 



and 



-e ut \ ( r \ _ ( ve ut p + e vt r 
i vt )\p)~\ e"*p 



(93) 



-e~ vt 



(94) 



F ={ 4> r J' B = -^{r' + rp + p' 

Take a simple transformation a as 



R = e vt p-e vto p, P = r-r, 

1 1 (95) 

R=-(r + r), p = __(e^+e 1 ' t ». 
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The Jacobian of a is 



a* 



1 
l 

2 






2 e 



-e^° ~ 





The inverse transformation is 



2 P + *< 



- e~ ut P, 



—e uto R-e- vto P, 
2 



and 



V- 







1 

2 



_1 

2 





\ 







consequently using (|95|) and (|97|) we derive 



dw 
~dt 



ve ut p + e ut p 
r 



-e ut r 



-\e ut P 
e—utR 



dw 
~dt 



\<r vt R- \i 

\& vt P + \z vt R 



Simple computation attains 



Set 



" 5 I =( ° 
^(i) = ^| _«*C°)^| =( -* Ut0R 

^ \t=t ^ww ^ \t=t I _ e -vt p 



W 



so we get a first order scheme for the system ()93[) 



r fc+i \ _ ( h£ 4+^ W rk 
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which is K(z, t)-symplectic. The transition matrix is denoted by A, then 



^ e --«+i o J \e*^* 
Similarly we get </>( 2 ) 

= ( ) . 

Direct computation also gives a second order symplectic scheme for the 
system (|93|) 

fc+l \ / 16-afe 8a \ / \ 

' ] — 16+a6 16+afc ] I fc I flf)^ 

n, , I \ 8a _-i/r 16-a6 »-i/T II „, I ' V iu0 J 
Pk+1 J \-T6+^o e T6+^fe e J \ Pk J 

where 

a = 2T-UT 2 , b = 2r + UT 2 . (106) 

The transition of the scheme (|105|) also preserves the symplecticity. While 
the direct Euler center difference scheme for the equations (|93|) is 



r k+l \ ( -t 2 +2ut+4 4r 

' I — I r 2 +2i/r+4 r 2 +2vT+i \ / 'fe 

-4t -t 2 -2z/t+4 

r 2 +2i/r+4 r 2 +2i/r+4 



(107) 



which is not a symplectic scheme. 
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